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Under many circumstances, the only important two-body interaction between atoms in ultracold 
dilute atomic vapors is the short-ranged isotropic s-wave collision. Recent studies have shown, how- 
ever, that situations may arise where the dipolar interaction between atomic magnetic or electric 
dipole moments can play a significant role. The long-range anisotropic nature of the dipolar inter- 
action greatly enriches the static and dynamic properties of ultracold atoms. In the case of dipolar 
spinor condensates, the interplay between the dipolar interaction and the spin exchange interaction 
may lead to nontrivial spin textures. Here we pay particular attention to the spin vortex state that 
is analogous to the magnetic vortex found in thin magnetic films. 
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I. INTRODUCTION 



Experimental realization of dilute atomic Bose-Einstein condensates (BECs) has revolutionalized the field of ul- 
tracold atomic physics [l|, Q, Q. For the first time, we have a macroscopic quantum object that is amenable not 
only to exquisite experimental control, but also to detailed microscopic theoretical description. Early experiments 
on atomic BECs were all carried out in magnetic traps, where the atomic spin is polarized by external magnetic 
fields, and hence the atomic spin degrees of freedom is frozen (see, however, Refs. [4], [5J)- The atoms, however, get 
their spin degrees of freedom back when they are trapped in off-resonant optical dipole traps [6], in which case all 
magnetic Zeeman sublevels of the ground state atom can be trapped. Such condensates are called spinor condensates. 
Collisions between atoms give rise to an effective spin exchange interaction [7j, [8[ , which is analogous to the exchange 
term in the theory of magnetism. The spin exchange interaction leads to interesting coherent spin- mixing dynamics in 
spinor condensates, a phenomenon that has been both theoretical studied [9|, LLC), LUU, LL2[ and experimentally observed 

Experimentally, spinor BECs have been realized in 23 Na and 87 Rb. The ground state of both these atoms possesses 
a small magnetic dipole moment of /i#/2, where \±b is the Bohr magneton. For typical condensate density (~ 10 14 
cm -3 ), this would yield a tiny magnetic dipolar interaction energy on the order of 0.1 nK per atom, which is a few 
orders of magnitude smaller than the total collisional interaction energy. Furthermore, it appears that the small 
dipolar strength would be completely overwhelmed by any finite temperature effect (typical temperatures in BEC 
experiments are ~ 10 nK). Consequently, it has long been thought that dipolar interactions in such systems can be 
safely ignored. Hence, experimental effort to achieve dipolar atomic condensate has been focused on other types of 

^il , atoms, notably 52 Cr [13, LL8[ which has a magnetic dipole moment of 6/^b in its ground state. 

C^> ■ Under a more careful inspection, however, the conclusion that dipolar interaction plays negligible role in alkali 

atoms becomes questionable. We realized a few years ago that dipolar interaction in spinor alkali condensates can 
play a more prominent role for the following reasons: 

1. Although the total collisional interaction strength is much larger than the dipolar interaction strength, the spin 
exchange interaction is not necessarily much stronger. Particularly, for / = 1 hyperfine manifold of 87 Rb, the 
dipolar energy can be as large as 10% of spin exchange energy, thus making a nontrivial contribution to the 
total spin- dependent energy. 

2. The spin-dependent interaction, although much weaker in magnitude than the spin-independent interaction, is 
a critical determinant of the magnetic properties of spinor condensates. 

3. The long-range and anisotropic nature of the dipolar interaction may further enhance its effects. 

4. A small finite temperature will not overwhelm the dipolar effects due to the fact that, in a condensate, the 
interaction effects enjoy a Bose-stimulation factor — the total number of condensed atoms. This same reasoning 
also explains the importance of the weak nuclear dipolar interaction in superfluid state of 3 He [19|. 

This motivated us to carry out a detailed investigation of the dipolar effects in spinor condensates. Our studies have 
confirmed that rich physical phenomena can indeed be induced by the dipolar interaction. 

In the following, we will first present the Hamiltonian that describes the system. We will then study the ground 
state properties. To this end, two approaches will be used. The first is the so-called single mode approximation — all 



spin components of the condensate are assumed to possess the same spatial wave function. Under this approximation, 
the Hamiltonian can be greatly simplified, which makes a full quantum mechanical study possible. The single mode 
approximation, however, assumes that the atomic spins are uniformly oriented in space. Its validity depends on the 
dipolar interaction strength, as well as other parameters such as the geometry of the trap potential. To go beyond the 
single mode approximation, we adopt a second approach — the mean-field calculation without any a priori assumption 
on the spatial wave functions. From this study we see that sufficiently large dipolar strength induces non-trivial spin 
textures in the ground state. The specific pattern of the spin texture is sensitive to the trap geometry. In particular, 
a pancaked-shaped trap favors the spin vortex state, analogous to the magnetic vortices found in magnetic thin films 
or disks. 

II. HAMILTONIAN OF A DIPOLAR SPINOR CONDENSATE 

We consider TV condensed spin / = 1 atoms trapped in an axially symmetric harmonic potential, 

y ext (r) = ±Muj 2 (x 2 + y 2 + AV) , (1) 

with A being the trap aspect ratio, and M the atomic mass. We have chosen the symmetry axis to be the quantization 
axis, z. The atoms interact with each other via both short-range collisions and long-range magnetic dipolar interaction. 
Under a uniform magnetic field B, the second quantized Hamiltonian of the system reads 

H = Ho + Hdd 
where Ho and Hdd represent the non-dipolar and dipolar part of the Hamiltonian, respectively, and are given by 



Ho = Jdrft a (r) 
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Hdd = 



2 j dr$,(r)$(r)^(r)^ a (r) + -| J dr$,(r)#»F a/3 • F a ,^(r)^(r), (2) 
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where F is the spin angular momentum matrices, e = (r — r')/|r — r'| is a unit vector, and ipa(r) the field operator 
for spin component (or Zeeman sublevel) a = 1,0, — 1. The collisional interaction parameters are [7, 8] 

47rh 2 (ao + 2a 2 ) 47rft 2 (a 2 - a ) 

Co= 3M ' C2= 3M ' 

where ap is the scattering length for two / = 1 atoms in the channel with total spin angular momentum F. Sym- 
metrization of many-body bosonic wave function dictates that only the symmetric spin channels F = and 2 are 
involved. The dipolar interaction parameter is 

with jjlq being the vacuum magnetic permeability, and gp the Lande g- factor. Finally, in Eqs. ([2]) and (|3]), and 
hereafter, it is assumed that the repeated indices are summed over. 

The first line of Eq. ([2j) represents the single-particle part of the Hamiltonian, while the second line results from the 
two-body contact interaction. The term proportional to Co is symmetric in the spin indices and represents the spin- 
independent contact interaction. The term proportional to c 2 , on the other hand, is spin-dependent and represents the 
short-range spin-exchange interaction. The sign of c 2 determines the nature of the spin-exchange coupling: negative 
c 2 represents ferromagnetic coupling, while positive c 2 represents antiferromagnetic coupling. The expression of Hdd 
in Eq. ([3]) follows from the dipolar interaction potential between two magnetic dipole moments /j, i = gF^B^i (i = 1,2) 
located at spatial points r and r 7 , respectively, 

T/ / /x Mo Mi • M 2 ~ 3 (Mi • e) (// 2 • e) 

^d(r, p ) = ^p • (4) 

The spin-exchange term and the dipolar term describe two types of spin-dependent interactions. It is the interplay 
and competition between these two terms that give rise to a rich variety of spin textures. 



III. GROUND STATE UNDER SINGLE MODE APPROXIMATION: SINGLE-DOMAIN STATE 

A. Hamiltonian under single mode approximation 

The total Hamiltonian of the system as represented by Eqs. ([2j) and (|3]) is quite complicated. There exists, however, 
a powerful method that can greatly simplify the problem. This is the so-called single mode approximation (SMA). 
More specifically, we assume that the field operators can be decomposed as 

2p a (r) =(j)(r)a a , (5) 

where </>(r) is a unit normalized spin-independent spatial wave function. Here we shall not worry about the specific 
expression of </>(r), which should be properly chosen to minimize the total energy. 
Inserting Eq. ([5]) into Eq. ([2|), we have 

Ho = J ' dr<f>*(r) (-^ + V ext (rfj <f>(r) N-g Ff i B B.t+^ J dr\4>(v)\ 4 N(N-1) + ^ J dv\4>(r)\\t 2 -2N) , (6) 

where TV = a^ a a a is the total particle number operator and L = d^ a F a pap is the total spin angular momentum operator. 
Following a similar procedure, we may obtain Hdd under the SMA: 



H dd = | J dvjdv' '^'J^p }l [(L 2 - 3(L • e) 2 ) - (27V - 3&t F a/3 • eF^ • eap)~ 
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(L + L 2 cos6> e sin<9 e e~^ e + h.c.) - ^(L-L z cos6> e sm<9 e e^ e + h.c.) - ^{L\ sin 2 6> e e~ 2 ^ e + h.c.) 
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where L± = L x ± zL y , # e and ip G are the polar and azimuthal angles of (r — r'), respectively. Obviously, this form of 
Hdd is still quite complicated. However, here we can take advantage of the spatial symmetry of the system to further 
simplify the dipolar part. As we have adopted an axially symmetric trapping potential, which is indeed the case in 
most experiments, it is natural to assume that the spatial wave function 0(r) possesses the same axial symmetry. 
Under this condition, it is not difficult to see that if we carry out the integral in polar coordinates, terms proportional 
to e ±imipe in Eq. ([7j) will not survive after integrating over the azimuthal angle cp e . We therefore have 



j dvjdv' |<K*W)| 2 1 ^l^ (~ t2 + 3l " + 3 ^ Q ~ *) , (8) 
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where no = %do is the number operator for spin-0 component. 

The total Hamiltonian under the SMA is obtained by combining Eqs. (|6|) and (|HJ). Since we are dealing with an 
isolated system, the total number of atoms is a constant. Therefore, we may neglect terms that only dependent on 
N. Finally, we have 0, HU 

Hsma = (4 - c' d )L 2 + 3c' d (L 2 z + n ) - gr^B • L , (9) 

where the two coefficients are defined as 

c' 2 = C jJdr\<p(r)\\ c' d = | j dvjdv* |0(iW)| 2 - ^gp ■ (10) 

It is worth pointing out the trap geometry-dependence of the effective dipolar coefficient c' d . If we choose <p(v) to 
be the single-particle ground state of the harmonic potential V^ x t, i.e., 




FIG. 1: The A dependence of c d with Cd — 1. 



where we have used the harmonic oscillator length y / h/(Mco) to be the units for length, the integral in c' d can be 
carried out exactly: 



CdVJ ' ^2A + l-3A tan " lv/A -" :1 



In Fig.HJ we plot c' d as a function of trap aspect ratio A. One can see that the sign of c' d depends on the trap geometry: 

{< , for a prolate trap (A < 1) 
= , for a spherical trap (A = 1) 
> , for a oblate trap (A > 1) 

This provides a convenient control knob one can use to change the properties of or even induce phase transition in 
the system. 

B. Ground state structure under the single mode approximation 
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FIG. 2: Phase diagram of spinor condensate in the absence of external magnetic fields, (a) Dipolar, (b) Non-dipolar. 

Hamiltonian ([9]) is reminiscent of the Hamiltonian that describes a quantum magnet. The ground state can be 
obtained by diagonalizing ([9]). In the absence of the external magnetic field, i.e., B = 0, the phase diagram can be 
plotted in the c' 2 — c' d parameter space as shown in Fig. [2](a). According to the nature of the ground state, we can 
divide the parameter space into three regions labelled as I, II and III. Region I represents a ferromagnetic phase with 
easy-plane anisotropy. The ground state wave function in this region can be written as |iV, 0), where we have used 
the standard angular momentum basis state |/,m) such that 



L 2 |/, m) = 1(1 + l)ft 2 |/, m) , L z \l, m) = raft |Z, ra) . 




FIG. 3: The ground state in Region III can be written as J2i 9i\l>ty- Here we plot the coefficients gi for c' d /c 2 — 0.2 (circle), 
0.8 (triangle), 1.4 (square), and 2.0 (diamond). We have taken TV = 100. 

Region II represents a ferromagnetic phase with easy-axis anisotropy. The ground state wave function in this region 
can be written as |7V, ±7V) with a two- fold degeneracy. In both Regions I and II, the atomic spins are aligned along 
the same direction, either in the transverse plane (Region I) or along the z-axis (Region II). Region III, on the other 
hand, represents roughly an ant i- ferromagnetic phase where atomic spins are entangled to form spin singlets. Here the 
ground state wave function has a slightly more complicated form: J^ 9i\l-,fy where the coefficients g\ has in general to 
be calculated numerically. Several examples are given in Fig. El The mixing of different / states is due to the no term 
in the Hamiltonian. In principle, the no term will have a similar effect in the other two regions I and II. However, its 
effect there is negligible for large particle numbers TV ^> 1. 

For comparison, we also present the corresponding phase diagram for the non-dipolar case (c' d = 0) [9] in Fig. [21(b). 
Under this situation, the Hamiltonian is simply H = c!^i? which possesses a full rotational symmetry in spin space. 
The ground state is determined by the sign of c 2 . It is an isotropic Heisenberg (anti-)ferromagnet if c' 2 < (c' 2 > 0). 
The effect of the dipolar interaction is therefore quite transparent: It breaks the rotation symmetry of the non-dipolar 
system and introduces magnetic anisotropy. 

Next we investigate the effect of a uniform external field. In particular, we are interested in the critical field strength 
at which the system is fully polarized by the external field. 

Longitudinal field — First consider a longitudinal field along the z-axis. It is easy to see that under this condition, 
L z is still a constant of motion since it commutes with the Hamiltonian ([9|). Its effect in Region II is quite obvious: 
Any longitudinal field will break the degeneracy of the ground state in the absence of the external field and polarize 
the spins along the field. Therefore the critical field strength here is infinitesimally small. For Region I, the new 
ground state should have the form |TV, m) where the value of m may be obtained by minimizing the energy 



E(m) = (N,m\HsMA\N,m), 



which yields 



m = 



Qf^bB 1 

c' + 2 



where [x] denotes the largest integer no larger than x. Critical field strength is reached when m = TV or 



B c = 



6c' d (N - 1/2) 
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For Rb condensate, this would correspond to a field strength on the order of 0.1 mG. The critical field strength for 
Region III can be obtain in a similar manner. Here the ground state has the form X^< mo ^l^ m o)- The critical field 
is given by 



B r 



(4TV-5)c' d + (2TV-l)c' 2 
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Transverse field — Next we consider a transverse field along, e.g., the x-axis. Its effect in Region I is similar to 
that of a longitudinal field in Region II: any transverse field would fully polarize the system in Region I. For Region 
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FIG. 4: The effect of transverse field in Region II. Classical picture of the spin orientation, 
degenerate when B < B c and fully polarized when B < B c . 



The ground state is doubly 



II, we can take a classical approach since here the total spin is macroscopically large (L « N ^> 1). Treat the spin as 
a classical magnetization vector with length L = N and polar angle $, we can write down the energy of the system 
as (neglecting the unimportant no-term) 

E W = (4 - c 'd) N<2 + 3c' d 7V 2 cos 2 tf - g F fi B BNs'mi9. 
Minimizing the energy with respect to $ yields the optimal value for d as 



Hence the critical field strength is 
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For B < B C1 the ground state is doubly degenerate as schematically shown in Fig.[H When B > B C1 the two degenerate 
states collapse into one and the system is fully polarized by the transverse field. We want to point out that, in a 
quantum mechanical treatment, the double degeneracy at < B < B c will be lifted by quantum fluctuations, resulting 
in quantum spin tunnelling [22]. Such quantum effects will be important when N is small. 

Finally, for Region III, the ground state can again only be obtained numerically as a superposition of different 
angular momentum states. The critical field can be shown to be 



B c 



^' 2 -c' d )N 
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C. Validity of single mode approximation 



Our discussion so far is constrained within the SMA. As can be seen, the SMA is a powerful approximation from 
which rich physics can be derived. Now we want to address the important question of the validity of this approximation. 

First we notice that the SMA implies that the spin orientation is position independent. To show this, we define the 
normalized spin vector as 

s(r) = S(r)/n(r) , 

where n(r) is the number density, and S(r) = ($£(r)F a 0^g(r)) is the total spin density. Under the SMA, we have 
n(r) = 7V|(/>(r)| 2 , and S(r) = n(r) (a^Fapap) = n(r)(L). Therefore, we have s = (L) which is spatially invariant. 
This situation would correspond to the so-called single- domain state of a nanomagnetic material [23]. The spatially 
uniform spin orientation minimizes the short-range spin-exchange interaction, but does not in general minimize the 
dipolar energy which is most effectively reduced for closure domain structures. Obviously the single-domain state 
is thus stable when the exchange interaction dominates. When this is not the case, the exchange interaction will 
be frustrated by the dipolar interaction, rendering the single-domain state intrinsically unstable to various other 
magnetically ordered states. Armed with this insight, we conclude that the SMA is valid only when the dipolar 
strength is sufficiently small and hence cannot cover the whole spectrum of interesting quantum spin phenomena. 

Our study confirms this conclusion. For large dipolar strength, the SMA is no longer valid and nontrivial spatial 
spin textures develop in the system. We shall now turn to this situation that goes beyond the SMA. 



IV. BEYOND SINGLE MODE APPROXIMATION: SPIN TEXTURES 

A. Mean-Field Equations 

To go beyond the SMA, we need to calculate the wave functions for each spin component explicitly. To this end, 
we adopt the standard mean-field approach. Specifically, the field operators ifj a in ([2j) and (|3]) are replaced by their 
corresponding expectation values ifj a = (^ a )- The Hamiltonian thus becomes the energy functional: 



Ebp a ,rj = J 



dr 
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(11) 



where n(r) = ^* (r)i/j a (r) and S(r) = ^*F a ^ p are the total number density and spin density, respectively. Min- 
imization of the energy functional E[ip a ,ip^\ requires 5E[i/j a ^^]/5ilj^ = 0, which yields a set of nonlinear coupled 
partial differential equations for the condensate wave function ^ a (r): 



ih 



dt 
dt 
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Vext - ^F^fB • Fq/3^/3 + c n 



*Po + [c 2 ^+ + c d P+]^i + [c 2 S_ + C d P_]^_i , (13) 



+ K x t(r) - ^fMfB • F-i/^ + c n - c 2 ^ - c d D 2 



^-1 + [C 2 5+ + C d X>+]^0 : 



(14) 



where S± = (5^ ± iS y )/V2 and the integral operator D^ is given by 

P,(r) = y dr' ^-^ [Si(r') - Se* S(r') • e] . 

The ground wave functions are obtained by solving these equations using the imaginary-time evolution method, where 
the term involving the integral operator T>i are dealt with using convolution theorem and fast Fourier transform. 

B. Phase Diagram 
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FIG. 5: Phase diagram of a dipolar spin-1 condensate with (a) ferromagnetic (ci < 0) and (b) ant i- ferromagnetic (c2 > 0) spin 
exchange coupling. Regions labelled by SMA-I, II, and III are where the SMA is shown to be valid. 

From this calculation, the validity of the SMA can be directly checked by comparing the wave functions i/j a . Phase 
diagrams such as shown in Fig. [5] can be obtained [24, 25]. Fig. 03(a) is the phase diagram for a ferromagnetic exchange 
coupling with c 2 < in the absence of external magnetic field. Here we fix the ratio c 2 /co to be —0.01 corresponding 
to the values of 87 Rb. The phase diagram is plotted in the parameter space spanned by the trap aspect ratio A and 
the ratio between the dipolar and the exchange interaction strength q = q/|c 2 |. A similar phase diagram for an 
anti- ferromagnetic exchange coupling (c 2 > 0) is plotted in Fig. 03(b). Here we use c 2 /co = 0.03 as in the case of 23 Na. 

From Fig. [5j we can make the following observation: 



1. For any given trap geometry, the SMA is only valid for sufficiently small dipolar strength, in full agreement with 
the qualitative argument we made in Sec. IIII CI 

2. In the region where the SMA is valid, according to the nature of the wave function which determines the 
local spin vector s, we can further divide the region into three subregions labelled as SMA-I, II and III. They 
correspond to Regions I, II and III shown in Fig. [2j 

3. The critical value of q where the SMA becomes invalid is sensitive to the trap aspect ratio: it decreases as A 
increases, i.e., when the trap becomes more and more pancake-shaped. This shows that the dipolar interaction 
plays a more prominent role in a pancake geometry, which is consistent with the observation in the study of 
magnetic materials — The dipolar interaction, normally weak enough to be ignored in bulk materials, plays an 
essential role in stabilizing long-range magnetic order in two-dimensional systems such as magnetic thin film 

M- 

4. In the region where the SMA is invalid, we can also divide the region into three subregions labelled C, S and P. 
We will discuss each of these regions in more detail below. 



C. Ground State Beyond Single Mode Approximation 
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FIG. 6: Ground state spin structure in (from left to right) C, S and P regions, (a) is a streamline plot of the local spins in a 
cigar-shaped condensate. Here the spins have a dominant ^-component and form a helical pattern. In (b), the arrows indicate 
the direction of the local spin vector, (c) is the spin structure in the z — plane for a pancake-shaped condensate. Here, the 
spins are planar and lie in the x — y plane. Spins at different layers of z have the same vortex structures. 



Now let us discuss in more detail the ground state structure of the spinor condensate in the region where the SMA 
becomes invalid. The phase diagram plotted in Fig. [5] shows that there are three regions labelled as C, S and P, which 
stands for cigar, spherical and pancake, indicating the geometry of the trap in which a particular type of spin texture 
is favored. 

Figure [6] illustrates the ground state spin structure in these three regions. In a pancake-shaped trap (Region P of 
Fig. [5]), spins are planar, i.e., lying in the transverse plane. Dipolar interaction forces the atomic spin to form a vortex 
pattern. We can write the wave function as ijj a (r) — \/ n a( r ) e lS(X ^ r \ Our calculations show that in this region, the 
wave function contains a nontrivial phase: 



6 a (r) =w a Lp + (j) a 



(15) 



where <p a is a constant phase shift whose absolute value is not important but satisfies (j)\ + 0_i — 20o = 0, cp is the 
azimuthal angle, and w a is then the "winding number" — circling once around the z-axis, the phase of the wave 
function ip a changes by w a x 2tt. The specific values for the winding numbers are 

(wi,w ,w-i) = (-1,0,1). 

Furthermore, we also have rti(r) = n_i(r). Therefore, we can see that the spin component 1 and —1 share the 
same density, but opposite winding number as illustrated in Fig. The spin vortex state is obviously analogous to 



®,(*,y,0) 




FIG. 7: (a) Density of each each spin component along the x-axis (solid lines). The total density is represented by the dashed 
line, (b) and (c) illustrate the phases Gi and G-i in the z — plane. 

the magnetic vortex state found in magnetic thin films |23|, |26j. However, there are important differences. In the 
central core region of a magnetic vortex, due to local magnetization conservation, the magnetic moment has to align 
perpendicular to the plane of the film in order to decrease the exchange energy at the center of the vortex. In a spinor 
condensate, however, there is no such constraint on local spin moments. Therefore local spin can simply vanish in 
the vortex core. Indeed, the core region is occupied by the spin-0 component while the other two spin components 
vanish. The spin vector can be expressed as 





= 2v / 2no^i 



S(r) = 



Such a spin texture is also referred to as a coreless skyrmion. Further, the spinor condensate is really a novel superfluid 
described by macroscopic wave functions which contain both amplitude and phase. In quantum mechanics, the spatial 
non-uniformity of the phase of a single-particle wave function is directly related to the velocity of the particle as 

v« = -V6.. 

We therefore conclude that, in the spin vortex state, spin components 1 and —1 circle around the z-axis in opposite 
directions, while the spin-0 component is stationary as wo = 0. This results in a situation with a net spin current but 
without a mass current. 

In a cigar-shaped trap (Region C of Fig. [5]), local atomic spins are predominantly aligned along the z-axis. In that 
way, most of the spins are aligned in a head-to-tail configuration which minimizes the dipolar energy. Here the phase 
of the wave functions can still be expressed as in Eq. ([T5]) but with a different set of winding numbers 



(w 1 ,w ,w- 1 ) = (0,1,2). 
Further the phase shift <fi a are no longer constants but functions of z. The spin vector in the C phase takes the form 



S(r) 



Asin((/? + S) 

-Acos((/? + 5) 

n\ — ri-i 



where A = ^2no(^/ni + y/n-i) , and S(z) = ipo(z) — ipi(z) is the spin twisting angle which is a monotonically increasing 
function of z with 5(z — 0) = 0. 

In between the P and C phases, is the S phase which occurs when the trapping potential is close to spherical 
(A « 1). A distinct feature of the S phase is that n a becomes non-axisymmetric, signalling the breakdown of the 
cylindrical symmetry of the spatial wave functions. A typical spin configuration in S phase is displayed in the middle 
plot in Fig. [6l One can clearly see that the spin structure does not possess the cylindrical symmetry and there is a 
180°-domain wall separates two spatial regions of spin. 



D. Spin Vortex in External Magnetic Field 



We have seen now that the dipolar interaction plays a more important role in pancake-shaped traps in which the 
ground state takes the form of a spin vortex. In this section we will focus on the spin vortex state and study how it 
responses to an external magnetic field. 
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FIG. 8: (a) The B-field dependence of magnetization M 2 for g = 0.5, N = 1.5 X 10 6 and A = 4 (solid line), 6 (dashed line), and 
8 (dash-dotted line). The inset shows the same thing except for q = 0.09 which corresponds to the real parameters for 87 Rb. 
(b) The A dependence of critical magnetic field for various (iV, q) values. 



Longitudinal field — First consider a uniform magnetic field along the z-axis. Here we define the magnetization as 

M z = 1 JdrS z (r) = 1 fdv [m(r) - n_i(r)] . 

We plot in Fig. [51(a) the magnetization as a function of the field strength B z . As expected, the magnitude of M z 
is a monotonically increasing function of B z [M z is negative here, because for the / = 1 hyperfine level of alkali 
atoms, gp = —1/2 < 0). However, there exists a critical field strength 5* at which M z suddenly jumps, signalling 
a first-order phase transition of the system. This critical field strength depends on the trap aspect ratio, the total 
number of atoms and the dipolar strength, as shown in Fig. Efb). 
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FIG. 9: The z (in units of ^/h/Moo) dependence of spin twisting angle 6(z) = <j>i(z) — (J)o(z)—7t/2 for q = 0.5 and N — 1.5 x 10 6 . 

To provide more insights into this phase transition, we study the wave functions ip a (r) — \/ n a( r ) e lQ(X< ^\ We find 
that the phases & a can still be written as in Eq. ([T5]) . The phase shifts <j) a always satisfy the condition <\)\ + </>_i — 20o — 
0. However, for B z < B*, (j) a are constants; while for B z > B*, <p a become functions of z as shown in Fig. [9l 
Furthermore, the winding numbers also change across the critical field and are given by 



(wi,wo,w-i) 



The local spin vector takes the form 



S(r) 



(-1,0,1), B<B* 
(-2,-1,0), B> B* 



n7 + y/n-i) cos((p — S) 
ri\ — n-\ 




with 5 = (j)\ — <j)Q — it/2. For B < B* , we have (5 = 0; while for B > £?*, S is a function of z as shown in Fig. [9l 

Transverse field — Now we consider the effect of a uniform transverse field along the x-axis. As the field strength 
B x increases, similar to what happens in nanomagnets [27J], the spin vortex starts to move away from the center and 
perpendicular to the applied field (along the ?/-axis), and eventually moves out of the cloud. This is illustrated in 
Fig.dOl 
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FIG. 10: Magnetic vortex under a transverse field for N = 1.5 x 10 6 , g = 0.5, A = 4, and, from left to right, the corresponding 
magnetic fields are B x = 0, 0.0143, 0.0286, and 0.0429 in units of mG. Shown are the spin structures in the z = plane. 

V. CONCLUSION 

To conclude, we have shown that dipolar interaction plays a crucial role in spinor condensates. We have provided 
here a detailed study of the dipolar spin-1 condensate both within and beyond the single mode approximation. In 
the SMA regime, the system is analogous to the single-domain state of a magnetic material with the spins uniformly 
oriented. For sufficiently large dipolar interaction strength, the SMA becomes invalid, and spin textures develop as 
a result of the interplay between the exchange and the dipolar interaction. We have found that a pancake-shaped 
condensate is particularly sensitive to the dipolar interaction. In such systems, the spin texture takes the form of a 
spin vortex which is analogous to the magnetic vortex observed in thin magnetic films. 

We have emphasized the connection between dipolar spinor condensates and ferromagnetic materials. We should 
also point out that there exist important differences between the two: While the latter represent a classically ordered 
system, the former are intrinsically quantum mechanical, as Bose condensates are macroscopic quantum objects 
described by quantum mechanical wavefunctions. Indeed we have seen that the spin vortex of a spinor condensate is 
intimately coupled with the topological charge of the condensate wavefunction: the spin components possess quantized 
vorticity or persistent currents. The fact that the spin textures and the motion of the atoms are intimately connected 
is a manifestation of an important property of the dipolar interaction, that is, the dipolar interaction provides a spin- 
orbit coupling. This can be more clearly seen if we notice that the dipole-dipole interaction potential (|4j) transforms 
as spherical tensors of rank 2 in both coordinate and spin space, and may be cast into the following form: 



Vdd(r,r') 




»S 2 



where Y2,m is a spherical harmonic of rank 2 and £2,™ its counterpart in spin space with components given by 



^2,o = \ - (Fi • F 2 — 3Fi z F 2z ) . 



^2,±i — =•=- (Fi z F 2 ± + Fi±F 2z ) . 
^2, ±2 = — - Fi±F 2 ± . 



The dipolar interaction is invariant under simultaneous rotations in coordinate and spin space and therefore conserves 
the total angular momentum (spin + orbital). However, it does not conserve separately the spin and orbital angular 
momentum. The spin-orbit coupling induced by dipolar interaction also creates the opportunity for the transfer 
of angular momentum between the spin and the orbital degrees of freedom under the constraint of total angular 
momentum conservation. This could lead to the Einstein-de Haas effect in the dynamical evolution of a spinor 
condensates [28, 29]. Spin-orbit coupling effect has also been predicted in ferromagnetic nanostructures [30]. 

Finally, we want to mention that the Berkeley group led by Stamper-Kurn have developed a wonderful technique 
capable of non-destrctive in situ imaging of local spin textures in atomic condensates [31]. Using this technique, they 
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have recently found tentative evidence of dipolar effects in spinor condensate of 87 Rb [32j- To summarize, dipolar 
spinor condensates represent a novel class of anisotropic superfluid as well as a new kind of magnetic material, many 
of its rich properties are just starting to be explored. 
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